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Under the assumption of constant failure rates it is pos- 
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elore te DUuLld a "relcability shorthand" which gives 
Diep URLeEted approach tO reliability computations for sys- 
+ems in the presence of complications like support by shared 
spares or changes in the failure rates of surviving compo- 
nents when other components fail. Th2> computational imple- 
mentation of the shorthand depends upon the convolution of 


strings of exponentially listributed random variables. 


This paper presents an algorithm for the numerical cornvo- 
lution cf exponentially distributed random variables. After 
reducing the system scenario to its shorthand format, one 
can use the programs that 2re given in the appendix ts 


obtain numerical values for the reliability of the systen. 
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Appendix A gives survival functions corresponding +o sev- 
eral reliability shorthand notations and a program writen in 


Fortran for computations from shorthand notations. 


Section III deals with the reliability of redundant sys- 
+ems under the assumption of constant component failure 


rates. 


Append#x® B consists o£ 2 program written in Fortran. The 
Paogream SUPPOEtS the approach of Section LIT. There isa 
crude Monte Carlo simulation program in Appendix D which is 


a simulation program parallel to the program in Appendix 8B. 


There is another program in Appendix C written in For- 
tran. This program uses the network approach to systens 


described in Section III. 


Appendix E summarizes the definitions used in this papér. 








This section introduces a general algorithm for con- 
puting the survival function of any convolution of 2xponen- 


tial random variables. 


In reliability, the «erm convolution refers to the summa- 
tion of independent random lives. In order to have simplic- 
meyem Specifying convolutions, the reliability shorthand 
introduces a special notaticn. In the following sections we 


Wieweuse this notation. 


Sune iVAuer ONETEON FOR A CONVOLUTIGN OF RANDOA 
TABLES 


hee a 
VAR 
ec F1(+t) and F2 (t) be the survival functions for the 

random variables T1 and T2 respectively. Let £1(t) and £2 (+) 

be he corresponding densities. Let F(t) be the survival 


hime tor f£O0r the random variable T, where T=T1+T2. 


Then the likelihood 2axpression for F(t) is 
L 
F(t)=F1(t) +f F2(t-s) £1(s) ds 
0 








[ieee wbedn= hand side Of €quatioa, FI(t) is the vrob- 
abllity hat component one completes the mission, 
E 
of F2ts) £1 (8) as is the probability tha* at some time 


s (OSsSSt) compenent one fails, component two takes its place 


mpomed=nces she system to fhe end of the mission duraticr ¢t. 


fo Order £O Illustrate consider som applications. 
1. Example 
Reliability Shorthand Notation : EXP {A,} +EXP {A2} 
SYSTEM: One component having one spare with a dissimilar 
maw bure races £& the active ccmponent fails, the 


Spare will replace it immediately. 


be seueMen life fOr the systemeis T-= Til + T2, The zeli- 
ability shorthand notation indicates that this system haés an 
exponential Life with failure rate A1 followed by an =xpo- 
nential life with failure rate A2. 
The survival function for the active component is 

= -A,t 

Py (T)= ¢e pee Ue 
The survival function for the spare is 

rot 


F fy el, (£20. 


The survival function for the system is 
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F(t) =A2/(A2- d4) git +)1/f (Ai ~A2) av 


Which is the well known result. 


Another way to establish this fornula is the use of the 
moment generating function. (Freund and Walpolef{ Ref.3]) 
2. Example 
Reliability Shorthand Notation : EXP{A}+F+EXP{ )} 
SYSTEM: One component having one identical spare. If the 
active component fails, the spare will replace it 


immediately. 


The formula that we derived abov2 becomes meaningisss, 


because the denominators become zero. If we proceed as 


before 
t 
F(t) = Pi(t) + fF2(t-s) f1(s) ds 
0 
= -At Cys = 
F(t) =e i ee 
0 
F(=) = (1+ y+) at Le eps ie 
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Tiiestesule Comes OL as expected 5 be the Bclang [{2,) } 


enavevel function. 


Be. RELIABILITY SHORTHAND NOTATION 


EXP {(A}+EXP {A} +....+EXP{A} (n identical 


nN 
mite =) (AE) S(i-1)! e7 =) 2). 
i=l 


The us2 of moment qenerating function gives the result 


immediately. 


Ce tHe RELIABILITY SHORTHAND NOTATION 
EXP {Az} wee | Ao} © sie eeeecee alee Ge {A} 


This is the expression for the convoluczion of n randen 


Variables where each has a distinct failure rate. 


By the approach used in Section 2.1, adding one exponen- 
mioweelrreras once, one can obtain th2 formula for the survi- 


Tater net. OT 


n 
ee ia, le mee een i, 20. 
i=l j4i J jAl 


| 
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This is also a well known formula and can be obtained 


from by a moment generating functior. 


D. THE RELIABILITY SHORTHAND NOTATION 


EXP {Ai} t+EXP{A2}+.. eaeeee50e eo et OKE TAD} (n, terms) 

Toe atop TORE { Nop tesccc a oo oe o EXD f 5} (n, t¢Frms) 

+ EXP (A, JH EXPLA, Ft eee eee coos e HEXP fr, } (n, +erms) 
K 

This is the convolution of } n. exponential random 
i=] 


variables where there are n, identical axponential randcn 


variables having the failure rate es 

The moment generating function technique is not useful in 
this situation, since there is a hug? amount of complexity 
mrvouwed., This section deals with this notation using the 


Convolu zion formula. 


EXP {A1} +EXP {A1} + EXP {A2} 


We know the survival function for EXP [)\,} +EXP{rA2} and 
also for EXP{A }+=EXP{A}. The next two subsections present 
fiereremen ways co reach the survival function for the shert- 


hand notation above. 


a3 





(a). 
Let T1,T2,T3 be random variables distributed as 


EXP {Ai}, EXP{A1}, EXP {A2} respectively. 


beet’ =T 1+ Zemernss random Variable has the Erlang {2,A } 


deste bution. 


The convolution formula for T=T1'+T3 is; 


t 

P(t) =F. (+) +| Fy, (t -s) za (cy ds ye se2 0) 

0 t 
Fy (t) = (1+)14) aie ee Pag oABas 

0 

= + =)ot . 
= (t) =e BU AyA2) 7 (ooo i)? vy AiA2oA (A 7A 1) C} 

“At 

els Rn (ene pce. 


(b). 


Meo iot2,73 be random varia gees distributed as 


EXP f{A1}, EXP{Ai}, EXP fA 2} respectively. 


Metem no =TO+T Ss. This random Fariabl> has the =eurvival 
mln G te. OT. 


7 At St 
F (=) = A» /{Ao7rA1) 2 UA Ome eS - Pat 


The convolution formulz for T=TI1+I2' is 


t 

Pr(t) = Fe (t) + [ F (t-s) f(s) ds, 220. 
2 ] 2 

0 
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as 


oS 
hoF (X27 1) ey On =>) ght | 
(ce eee 2) ds 


9) 
“Ai 
= 


ie 
((A2 2-22 3A2) S(A2 “Ar ) 2 + Ad ZF (A27A1) 3 
e2t 


Gee On ane fp 20 
Subsections (a) and (by) 


illustrate that the convolution 
formula gives a unique resul+, regariless of the way of 
Snoosing the prior random life. 
Pe 


EXP {A\,} +EXP ()\, J+ EXP (AX .} + EXP Q, } 
bev TA, 


T2, T3, T4 be random variables distributed 


Pema}, EXE {A}, EXP {A2}, EXP{,,} respectively. 
Then from the derivation 


Sect Onme2 se Liat ie tee 38S 
random veriable having the survival function 
Fry ({) = (a ta I Se ES) 
;' Ii £2 
where 4a11=(\¥- 2),)\2) 4 (do 7 Ar) f 


0 ALF KOT ) » 


The convelution gormula 


fom L717 bt end 
* 


is 
F = Pre fe F = £4(s Fer > 
Fy (t) SGA + JE ee Ss) xi jds 0 


a -X 
Pr (t) =(aiitaizt)e 


Teas ete (ale ( Hh, Ar &2tA 12128) 2 se 
5 0 
+ i> Se Sls 
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The result from the above is 
F(t) = (attaray Ot (age aye) EME, 220, 
ere ns aoa 7 Oo Ay) = eto Aaa <7 > -A,) es 
Senne Ss oy (lero) ge A < to UA, As)? oe 


It 1s important to nore that the number of exponential 
terms equals the number of dissimilarc failure rates and each 
exponential term has a polynomial co2fficient with the 
degree of the polynomial equal to the number of identical 


random variables having the corresponding failure rate. 


The next section deals with the convolution of exponen- 
fiaiy d:StrLpbuted random variables ising the fact iilus- 


rated above. 


Peomtesubsectien 2, we Can infe= the fo=n of the sur- 
vival function which we sought at th2 beginning of 


Seats on, C. 


SHORTHAND NOTATION : 


EXP fA 13 tEXP fAy } + eee coe EXP { A; } (ay) 22615) 
+ EXP {As} FEXP [Aa } + eee eee EXP (Ad} (mn, terms) 
+ EXP {AL} + EXP {AL} teoee ee EXPLA, | Ghie SENSENS}) 
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SURYIVAL FUNCTION ; 


F(t) = Ai(%) oe, Rea) unis (+) at Be) 
ice 
where A. (t)=a11tai2 tha tet .... aon i fe 
nol 
A2(7)=a21+4 get a 290 = tiie ae ; 
2 
n,-| 
k 
a. Examples 
SHORTHAND NOTATION : 
EXP {A,} + EXP fro } 
+ EXP {), } 
+ EXP{\,}+EXP{A;,}+EXP Qj } 
pumveval FUNCTION. ¢: 
F(t) =(a11+a12 ty eit neato grok sas, wel ay) c+a33,¢2) eo?" +20. 


7 





a set igortthm © eo Compuce “he Cock ticients 
The algorithm represent2d below develeps che 
survival function by adding one random variable in each run. 
As an example , in order t> compute the survival function 
for the convolution of ten exponentially distributed randon 


variables, the algorithm is supposed to run ten times. 


The notation used in the algerithn is: 


K number of dissimilar failure rates 
A; ith type failure rate 
d, Maa lume pete for the currently entering life 
2 
airy k*h coefficient on the jth polynomial 
n curren+ number of identical lives having 
the ith failure rats 
nn. number of randem variables having ith 


failure rate. 


Ii Se leereedr. i Bl 2 3, =0 » VK where j=1,2,..K , KEV, 250-62; 


EO lice di, pee Whe Sete — 0 2. aK 


Where = 0770.1... K 


J 
ee 
= 

= i 
4 


Riemer Sseaerun is °: ny=al, a1, =1 
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iigOraenm — 


1. Update the coefficients : 


So a , FOr K=2,3,--,0, 
2 
a. ppA a. /({m-1) , where =n, =) 0 Or j=9,1,-,2, - 4 
e e em-| e e 
2. Update the coefficient a 
n. & 
i 
: J 
a; az qt) (ahi Ged) +) (F-1) 14 ag s7(A2-Ad 
"6 'e! if, ut "e j=2 a tJ Ve 
n,>0 n,>l 
Semucdace the other cost ficients: ay Vay k 
where ~#2 ,n. #0 EDC a; 2 ene K=1,2,.,2, 
@., =a. A. f(A. A.) , Wi where ifi. 
in. ins i, ee 
= = = j S¢ 3 #7 i Mice 
2a ant, mas yi 7 OA; A) Pevlewnere sf: and ni2) 


: 
Or m=n,-J, J=1,2,--,n,-1 
Se. Example: 
Reliability Shorthand Notation 
EXP (A, }+ EXP {,, }+EXP {,,} + EXP h, } 


meter iet2- 13 Te be random variables distributed as 


Pump eiye eke ih,j}, EXP {ho}, EXP{,, } respectively. 
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VemvecumGdelikouwtOeGer.V= a f£O=mula £55 =he survival func- 


eeon or the random Variabls T=T1+T2+°3+TU. 


The use of algorit hn: 


1St Roun: 


At the end of this run we have only one random variable, 


meen =s Gistributed EXP{ A,}. 


—_ oe = 
2 ass then @1, 7 hy 
a re — 


iemeeee end Of the 2. Fun we have the stt®rviv2il functisa 
ror the T'=T1+T2, where T1 and T2 are identically distri- 


buted as EXP{\,}. The survival function is 


F (=) = (a, y+a.%) eee 


Which is ERLANG (2, Ai}. 


3% ROO: 
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Soe ray ian (kgs As ) ta UE 5s 7 (Ag = Ag) = 
=TA1 7 (Ara 22) + AtA2/ (Al 7A2) 2 


=177 (A17A2)? 


Q . =@12\2/ (22 7A1) then 412 FAAe Og ~ ry) 


211 = {Ao Qi:74)/7 0. se Shen te oe eT eo as A 


Note thet, here the coefficient aiz2 is the updated one. 


At the erd of the 3. run we have the survival functicn 
for the random variable T"=T1+T2+T3, where T1,7T2,T3 are as 
defined before. 


Ea -\,t -Aot 
Fait) = (a1 # 2b) € taoe 9, &20 


Where a,= 029-2 AD7 Og WALI Fe ALFRAAS (AQ AWAD) # 
aor = A177 Oats 
Note thet the coefficients are identical to the result of 


Subsecticn (>). 


ut) pon: 
= 726 n, =2, n, =2 
Qo2 =A2%21 then aa eens) * 


Soon yh 7 (Ay TAD) eiee/ (= AD? 


then aime ee ss 1? Wee) 





ae aoc vd Pim) ~ oven Q@12 =Arrj2*7 Cd27A, ) F 


enue eee pe iy) Nee Nae) 7 (2), ) 
At the end of the &. cun, 


we have the survival function 


meee the Variable T=T1+T2+T3+T4. 
F(t) = (arta?) sty (a, ta, ty H2t , +20 


Where a1) =(Q2 373A 7)2 =) 70027-0107 6 8842 ie? Oona f 


a@o1 =) Pe a he re a) , 84, Ni Aa SCOR) © 


Note that the coefficiants are identical to the result of 


Subsection 2. 


Ze 





This section deals with a system whose components haves 


constan* failure rates. 


Having the reliability network for a system in which each 
component has an exponential life and knowing the propabili- 
ties for failures (discussed in Appendix £2) makes it easy 


to describe the system's life. 


In order to make the idea clear, we will go *hrough some 


examples. 


Ae. EXAMPLE 
a 4 A2 


Vy 


2 spares (A1) 1 soars (A3) 
SYSTEM : 2 components in series with spares. 


DESCRIPTION : System has 2 componants in series. Compo- 
nent +t has a life distributed as EXP{A1} and thare are 2 


Glen c2celesparcs. Componens 2 has a life distributed as 
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BXP{A2} and there is one non-identical spare, whese life is 


Gesteibuted as EXP{)\ 3}. 


las 7 
Be eas 
3 AL /(A1 +2) 
as smal EX? (st 
ete ths ) 4 oS Oy *A2) LD ORAL ee, 
EXP {}, +2} \ ae: +, } + 
Cea Xo) Mee ys aa 


EXP Qi tag + ee cs} 


Seen (og bo) 97 Or +3) 


asia ee. ty, 3+ 
10 A347 (dA, +3) 


ZERO 
EXPplanation for the derivation of system's life: 


At the beginning the system has an EXP{\; tA.} Life. The 
failure of the system is by the failure of component 1 with 
Eee Oosoiietty Of A1/(Ar+A2) Of by the failure of component 2 


Meer oma ds lity Of A2 / (jy tA); 
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[ioe nee lite Eigqure, = number “@deistes the even] "Lailire 
of component 1" and number 2 denotes the event "failure of 


eemponent 2". 


If event 1 occurs, component 1 is replaced by one of ths 
Spares and system again functions with a life distributed as 
EXP {A; +A2}, Since the exponential distribution has the mémo- 


ryless property and componen* 2 still has the same failure 


T&F the event 2 occurs, component 2 is replaced bv its 


Spare and the system has a life distributed as EXP{hk +\3}. 


Mine Lumpers Of the lifes figure cotrespend to the transi- 


(D 


Geonms that Can occur. Th probability on each arc shows the 
conditional probability of the transition. As an 2xample 


eememmcam OCCULD Wit hwprobabil=tyyoft A: /0,,+),) pe gtven chat 


eyent 1 has occured before. 


The distribution ZERO defined by Esary[Reft.1j] and 
Repicky {Ref.2] (also defined in Appendix £3) enters when 


life is exhausted. 


rh 


ine 


For convenience of description, it is helpful to de 


the concept of path used in this paper. Path denctes the 
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Sequence of events in the system's life from the starting 
point to the point where the system is not functioning. 


Examples; 


Events 1 and 3 are a path, which denotes a sequence of 
lives for the system. In this cass, the system has 3 e2xpo- 
nentially distributed lives EXE Ary eke (ee Pee he). S 


Reo che probabilicy of Dy AN IONE aN oe 


Events 1, 4 and 8 form another path, which describes 2 
sequence of lives for the system. In this path, the system 
Meroe ss EXP (NG #7} += XP{), #)> FHEXPE LA, + A } +ZERO. The prob 
De veOrom= nes Dath is 


rae +o ih. AON + rallrs SQ. +h )- 


PiemeonO da sceri button Contra bU ss zero additional life to 
she sSyYS-em, SO we Can omit it. Nevertheless, we can not omit 
PeseoraQbanpll2cy in the celcuiation of path probabalety, 


Since the event numbered 8 has a probability of occurring. 
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Spemeeer SURV aA FUNCTION FOR THE SYSTEM 


Paria s Seecclon We Wail deal with an appreach t9 obtain- 


ing the system's survival function. 


SYSTEM: Two components in parallel. 


r4 








hy Z4A1 F Ad) EXP {i,} 


EXP {Ai1+ 2} + 


\ 


ho 701+ Ad) SXP {ii} 





Des eretON: Component 1 has EXP{\,} Life and component 2 


has SXP{A2} life. 


There are two failure events that can be allowed. The 
mie st as chat component 1 faiis at some time * and component 


2 carries the system for the rest of the time. The cther ona 
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2S that component 2 fails at some time during the mission 


aeascompenent 1 camries the system for the rest of the time. 


There are more events such as no failures during the mission 


duration, which are taken care of by the representation. 


Now we have two paths 


No.of Path Weighs Lee 


cok 


oe he ON PAPAS Ap PP exe 
Zz nS a a EXP ie +), }+EXP Get 


Pe steer eneaseStem’Ss time to failure ard le+ 1T1,T2,T3 be 
random variables exponentially distributed with the failure 


mares Nyt Xoekrs Ay TeSpectively. 


Then 


[1 + 72 with probability Ps 


Ti + T3 Wie Ons OU ety 95. 


BPiemsucveval fEuncticn can b2 written as 


yeep P(t) +p F(t) » 220 


Meso ye, (A itAa)s Po = \o7(y FAn} and F, (4), F, (t) denote 
Siomaweyry a  cunecticens for the shorthand notations 


EXP {A, +A. } t+ EXP{A.}, EXPLA, +} FEXP (A, } respectively. 


28 








The survival function for the convolution of two expcenen- 


tially distributed random variables with dissimilar failure 
rates is 


Ast 


MEE) 6 Dey take) G2 e200. 


If we do the substitutions for F, (t) and Fo (t) as 


Ay =Ai+h2o,A2 = Az2 anddAr =Aritr2 ,A2 = Ay Tespectively, 


miele o =cOm®es 


as a ote a + 

ey =e (= 1) Ane A2lt oa, +A2) SAi ee Pe al) 
and Rate) becomes 

os ~(Ai+A2)t -rot 

FoR) FA 7 (m2) 2 Tee We Xs © , +20. 


Then 


P(t)=Ar7 (Ai tA2) Fy (2) FAo7 Ay FAL) F (2) 
F(tjse! +e4} afritha)t 5 Ee 
The result gives the survival function «hat is well known 
EOE thas SyYSten. 


Ancther example is the 2 out of 3 system. 





RELIABILITY NETWORK: 


fijemicnownh SUrVival function is 


7 —(ArtAg)t -(AyeAg)t -(Artrg)t _ -(Aitr2+A3)t 
F(t)=e +¢< +2 -2=2 ; <20. 


_ 


it fs 


BX P{\, +); } 
Ai FZ (Ay FAQ FA3) 
EXP {Ar trot Ag} +—A2Z7 Or tr2t dg ) EXPLA, +33 
As7 (Aa tA2 +A3) 
arcs casi 
EXP {Ar trA2tA2 } FEXP {Ao+)3 } with probability \V7 (Argh) 
Meme CPt tls } FEXP EX, +A} Weeh probability ),/ (\,4,%,) 


Bee ato P+ EXP {i +45 } with orobability 7 (\,**,) 
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T= SUTVival function is 


F(t)=p; Fy (t)+po% (£) +p, Fa (t) » ¢20. 


where 
= -(AitA2tr3)t -(AztrA3)t 
Fi (=) = QotAd/ (-Xx) @ + Ar Hz #A3) 7A, 2 
= -(Aytr2tA3)t -(AtA3) t 
Pee ome stig eT egy A, oO 
= -(ArthotA3)t -(A1tA2)t 
FCY=QHI/ (ay) © FOAAIASe 


If we do the necessary cancellations, we can det, 


gc (AztAa)t p-Outas)t (Ar +A2)t (ArtAatA3) t 


F(+)= +e eu , t>d. 


as desired. 


2. General Procedure 


Having the algorithm presentei in Section II, we can 


treat more complicated systems under similar assumptions. 
The procedure is 
i. Set up the reliability network for the systan. 
memeneceording © this nStwork, set up the systen lige. 


Meewecrmg tne proper Feliability shorthand formulas for 
sierecelaced Convolutions Of exponential lives, set up the 


Survival function. 
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Ife SUA BY 


The reliability shorthand is an 2asy way to describe a 
Swomemecw lee Dut it as ditficult to 2mplemere ccomouta ioc- 
Gaye Saurc] there is considerable comolexaty in handling con- 


VOLMe ONS. 


The algorithm presentei in this p2per gives some relief 
Peomethnis G1fficulity. However, the accuracy in sbtained from 


Misco adlqg@micht is very much related to the differences in 


MmeGene= asp=cc in the al gorithm 2s that distributions are 
convolved one at a time and this requires very 2ccurate run- 


Ning conditions in the casa of a complicated systen. 


Tt is believed that it 1s possibls to derive another 
eugrerithm whach is more psyerful than the one introduced 
hers, Instead of adding soe distribution at a time, one can 


Gaeomeaaq Several at a tine. 
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A- 1 
ioc see  OnMCOntatns =<he SuvvVival functions for several 


shorthand notations which were derived by «he use of thea 


approach described in Section II. 


Aa latesmem@ehand Notation : SXP{ A}. 


SuevaVa ! Function: F(t)= es Mo cae 


Remipeeonmsmorthand NOohat: On: EXP{A,}+ EXP{)\5} 


Suse veva bk FuAction: 


F(t) = A2/(A2-A1) 54 t AGS (A. wAz2 ) shot peso. 


mee SnO=thand Notatleor + EXP{A}+EXP{ A}. 


SEV: Val Function: 


F(tpe (t+ retype ** | +20. 


erie onostnand NOcd ton 7: EXP (),} +ExXP(,,)+=ExP 4h, } 


Silay va0 Functions: 


F(+)=ai1 SU ea Be 


+35) Se 0. 
Premera n.7 (4, %) s7a,)* a, =AiA3/ (A, wA2) (A3 Az) 


a3 





Bay eA (Aya A3) (A2- X32) = 
Pein S oe?) EXP (A PFENP{ A} * EXPLA } 
SMevival Function: 
Beye 17202) a tO. 
A. 1.3.3 EXP{\1 } +EXP{A2} + EXPf{i2 } 
Surveve | Function > 
(ce ae tye es) 
Where a= A2*/ (dA27A1 Ae B= (Ay 27 2 A2) 7 (A 7 A2) ? 
Bo AyAH (A172) = 
Aw 1.4.1 EXP{A1} +EXPQM([2} + EXP[LAs} +EXP Ly } 
Sumvava Le Punction: 


ee ae! taper? 4a, 63 toa, ot Pie ea Oe 


where ai, = Il Aj/ Il (Ap Aq) Vota tae 
j#i il 
Aw1.4.2 SXPLA }+EXP({ AP+ EXPLA J + EXPL A} 


Survival Functior: 
P(t) = (14924172 4 2t24176 1343)8*" | t20. 
A. 10.4.3 EXP{\1} +2XP {ro} + EXP(A2} + EK? {ro } 


Sueveuvalo eunce on: 
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= -Axt -r2 
Mec = tf {a,+e,tta t<je , = 20. 


Wineme 4@.= Ns 2742-2 4 =e a 
aoo= Manele 5) 7 (Xk ae, 2 ’ ade Orie Wea oi Soe? A 
A. 1.4.4 EXPO, JFEXP QA, } + EXP{A2} + EXP Lr } 
Sieve Val Furction: 
= -A2t -Ast 
F(t)=(a,,ta,t)e gr ((slph ze At) SY, eralie 
where a = Cr29-SAFAL) SZ O27 AL) Fe 8 ALAEZ (Xn 7 A) 2 
an 0, J Sa, 7 Ure oo) ie Ao FAnki 7 (Ay 7h ) 2@- 
Ae1.4.5 EXP(A] } + EXP {A2} + EXPD{A3} + EXO [dg } 
Survival Fuaction: 
ae -Ait -A2t -r3t 
F(=)=a e@€ +a @ +(8 #2 t)2 ee Oe 
ll sa | Sct 32 


where a = A2\37(,o--1) (smh ) 26 ay = AZZ (Ale) (Ag "he ) 2 


Qs pho OCArrAg) Cag 7 Ag) tAaAaAgl IZ OA 72) (Ar 7A3) 


SAO a OE Cn) a 


a = W2h3/(,- Ag) X27 A3) = 
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A=Z 


This section introduces a Fortran program using the 


algorithm described in Section II. 


Cc 
C THIS IS A PROGRAM TO COMPUTE THE RELIABILITY OF A SYSTEM 
C WHICH HAS THE RELIABILITY SHORTHAND NOTATION 
C EXP (L1)#. .+EXP (LN) 
C WHERE THERE IS NO RESTRICTION FOR THE FAILURE RATES. 
c VARIABLES: 
C A(I,J)  : I.TH TYPE FALLORE RATE, J.TH COBFFICIENT 
C ON THE POLINDM. 
or M(t) > AMOUNT OF LIVES HAVING THE I.TH TYPE 
C FAILURE RATE. 
C NINIT(I) : AUXILAURY ARRAY FOR THE NI(I). INPUT 
C FOR THE PROGRAM. 
¢ 2d) : THE AREAY FOR THE FAILURE RATES 
REAL L(2C),A 10 
: INTEG? Nehtd 43.6). hr (20) 
C GET INPUT 
READ (5,499) T 
WRITE [6 498) T 
READ ( 500) K 
WRIT ge 1500) K 
READ ( 2301) (NINIT (I) ,L(I) ,I=1, XK) 
WRITE (6,1501) (NINIT(1) ,L(1) ,L=21,X) 
C 
C COMPUTE COEFFICIENTS ONE ALL 
E DISSIMILAR 
KK=1 
2 A(RK, 1) =1. 
NI (KK) =1 
hee 
4 TE (JJ +20; KK) 39 TO 
A (RK, 1) =A (K 71) *L (3d) / (L(Jd) -L (KX) ) 
IF(3d.=Q.k) GO TO 9 
JJ=3d+1 
GO TO 4 
9 IF(KK.EQ.K) GO TO 8 
KK=KK+1 
GO TO 2 
~ 8 CONTINUE 
C BEGIN T) ADD ONE AT A TIME 
C 
IP=1 
32 EE (NENIZ (TZ). EQ.NI(IE)) GO T9 93 
‘ NI (LE) =Nz (Iz) + 
€ UPDATE If 
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This section contains a program described in Section III 


mOomecOMpUuve mene reliability of 2 system. 





C 
eC 
: THIS PROGRAM CALCULATES THE RELIABILITY FUNCTION: 
C 
Cc F(T) = 
Cue? Ween (=P XP 4nd, 1) 4EXP(L1,2) +....00.tEXP (11,81) )} 
@ie PQ P20) (=EXP (LZ, 1) tEXPt(L2, 2) +..0. 00 etEXP(L2, M2 
c * PouetPeerait) (=fXP (19,1) +EXP (L3,2)+..00000t2XP (LI, M3 
+ @ ® e @ @ @ @ @ e e a @ a) e e og @ @ e 
cS @ @ @ @ e e = e £ Ss e e @ @e @ e @ @Q =) @ 
C + @ @ @ @ e e @ @ @ @e e @ @ e @ e @ 
: + PI * FI(T) (SEXP (LI,1) ¢2XP (LI, 2) +....... EXP (LI,M1) ) 
+ @ @ @ @ @ @ @ @ @ @ @ e @e @ @ . @ 
C ae @ 2 @e @ @ @ a @ @ @ @ @ e @e @ a @ . J a 
C 2 e @ e s @e e @ @ e @ @ @ @ @ Q @ @ @ @ 
C + @ @ e td @ @ @e e @ @ e e @ e & @ a 
C + PN »* FN (T) (SEXP (LN,1) +EXP (LN, 2) +.0e000e+tEXP (LN,MN)) 
C L > ARRAY FOR ALL LAMDAS. 
C P > ARRAY FOR ALL PROBABILITIES. 
C MI () - NUMBER OF FXD IN EACH ROW. 
C IK - TOTAL NUMBER OF ROWS. 
C PROBA : PROBABILITY OF SYSTSY SURVIVAL AT TIME T. 
Cc TIME > TIME 
C A(I,K) : TH2 COEFFICIENT FOR [THE CURRENT PATH, 
C I.TH TYPE FAILURE 
c RATE AND K.TH COEFFICIENT ON THE THIS POLINOMIAL 
REAL P(50) 
CALL READ1(IK,P,T) 
SUM=0. 
DO 1 I=1,IK 
1 SUM=SU M+ (T) 
IF (ABS (SUM-1.).GT.1.9E-5) GO TD 199 
PROBA=0.0 
DO 2 I=1,IK 
WRITE (6, 1009) sages 0) (00) 
CALL ONECON(PRO.T 
2 D ROBA=PROBA+ PRO*P ( I) 
WRITE (6,601) T, PRO 
STOP 
199 WEITE (6,198) SUM 
601 FORMAT ('1", 10X, "PROBABILITY ( TIME > ',F10.5,")= ', 
1009 PORMAT 10°, 10X,12,'.' , 3X ces ee 4p F10.7) 
198 PORMAT('1',10X,'THE SGM OF THE PROBABILITIES IS NOT 
* EQUAL TO 1.0°. /10X,'SUM = ',?10.7) 
e 
A GET INPUT 
C 
SUBROUTINE READ1(IK,P,T) 
REAL P (50), 
READ (5,499) T 
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INTRODUCTION 


This section consists 9f a computer program to compute 


SyoceMeccilabilz “ies as described in Sect®on Ir. 
again for simplicity, we will go zhruogh an example. 


The structure of the system is 


spare (A1) 


Ag 





e 2 
Components 1, 2 and 3 have lives axponentially distributed 
with failure rates ri sro 7A; respectively. Also we have? a 


Spare £Or component 1 which is identical to component 1. 


Using the shorthand approach, the system life would be 


detezmined as follows 
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Ay Z AitAztA3) —EXP fAAA3} 
EXP fth As} t—A2 7 AithatAJ—EXP fy 43} 
AiZ O,%, *,) A3/ (i, th, ths) —-ZERO 


cr. 


EXP {\,+A,+,} sill +),) -ZERO 


of 0, +2,+2,) A / Oy ty)—EXP 043) 


mis EXP {A+} + 
A3/ (\, +A.) —-2 ERO 


Sceme definitions are necessary before describing +he pro- 
gran 
NODE: Fach nod2 represents an exponentially distrubuted 
Eandomivecieble with a certain failure rate. The number for 
any nede can be chosen arbitrarily but can not be used nore 
than once. 
AEC : Bae wcre (Orig natss et a mnod2 and leads to eno=her 
Medes Only Ohe 2rC Caf tal Binate at a given node. Arc nun- 
Pemomeew be Chosen arbitrarily. Ifo. is the total number of 


Bet@eser 2 System iife than there will be n-1 arcs in <his 


system. 


tt 
iD 
4 
=| 


With these definitions, we cin reprasen+ a system 1: 


the following way 








BACK BOANTER LIST (IPB): Bech eped= has a back pointer. A 
back pointer is an arc number which shows which arc connects 
the node to the tree. Phe ne meOOTa ces 2S Zero, then the 


related node is the root of the tree. 


NODE NO BACK PIINTER IPB(T) 
1 0 
2 1 
3 ? 
u 3 
5 4 
6 5 
7 6 
g 7 
9 g 


(Here node 1 is the root of the tres.) 


OWE ist: As we mentioned, 2ach node represents an 
exponential lifetime, for simplicity we can use some integer 


code numbers for each failure rate. 
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a Le _ a= 


0 (aa Serabueion ZERO) 
1 yo eee 
Z | + Le 
3 Sh ae 


Bememenmeyectcrs Code number 0 is not arbitrary and 0 can be 
WeeaqmsOn fhe ZERO distribution. The sthers can be picked out 
Gaia trar ily . 


AR 


CObE>L Ist: This list is similar to ths node code list. 


Arc code numbers represent probabilities. 


Code No. Probabslity-PAa (TI) 

1 Ai 7 (Ay tA2*0 3) 

2 N27 (A1 FrA24*A3 ) 

3 h3 7 (A 1 FA2tA3) 

4 Ay 7 (1 tA 3) 

5 X37 (A *A3) 
ReGwoumeiN LIST : BPaewealcenhas am -ob2g2n node and @ temni- 
nal node. In the pregram we need to use only the origin 
ese . 
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Arc No. Origin Node £0 (i) 
1 1 
2 1 
3 1 
uy Z 
5 2 
6 2 
7 3 
8 3 


Peo bebo unn NODE “E@oT {LP}: This list indicates the nodes on 


the end of each path. Th2re is no necessary order in the 


mest (eee DEPTH LIST (19D): This list indicates the number 


OL vanesce trom las= nede <0 the root of the =ree,. 


i &- 2. Node List LP (Z) Denes De Pe HLS Seer) 
1 5 2 
2 6 2 
3 ¥] 2 
4 8 2 
5 9 2 
6 4 1 
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a 6 
(47-20. 0) GO T) 3 
} =L(LCO DE(I) ) 


NE.LCODE(J3J)) GO TO 4 


CONTROLS THE 
oo SUNG I EON OF PUR CURRENT PATA 


ih RO, TY 
L, PRO, TI 


pees ONE CONPUTES THE COSPFICZENTS 
Pomened IS DISSIMILAR 


SUBROUTINE ONEDIS (Ky A, NI, L) 
REAL A(11, e1Shs L (11) 
INTEGER NI 
a (KK, 1) <1 
NT KK =1 
JI=1 
TE (J5,2Q.KK) GO TO 6 
A (KK, WHA (KK, 1) AL (SS) / (L (TI) ~L (KR) ) 
TF (JJ.EQ.K) GO TO 9 
JI=TI +1 
GO TO 4 
IF (KK.EQ.K ) GO TO 3 
KK =KK+1 
GO TO 2 
RETURN 
END 
BEGIN TD) ADD ONE AT A TIME 
SUBROUTINE ONEATA(K,A,NI, L,NINIT) 
REAL L(V) A111 } 
INTEGER NINI (14), NI (11) 
EF (NINIT (TE) © EQ-NI (TE)) Gd T9 99 
NI (IE) =NI (Z5) + 


UPDATE I 
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ANA AAAQ AANA AN 


ANAND 


AMA 


AAD 


GET THE TIME 
READ (5,505) TIME 

GET THE NUMBER OF ROWS 
READ (5,501) N 

GET ALL P'S 


READ (5,504) (P(I),I=1,N) 
GET ALL LAMDAS IN THE ORDER OF BON 
MIDO=1 
MIUP=0 
DOmOOpJ=1.)N 
re (J. NE. 4) MIDO=MIUP+ 1 
MT (J) =MIDO 
READ (5,520) MI(J) 
MIUP=MIUP+MI (J) 
READ (54 50 3) L(1I) ,I=MIDO, MIUP) 
100 CONTINOGE 
RETURN 
501 FORMAT SE, 75) 
503 FORMAT(5F10. 
504 FPORMATI5F10.7 
505 FORMAT(5X,F10.3) 
520 FORMAT (5xX,15) 
END 
SUBROUTINE CONTRO(L,N,TIME,MI,MT 
REAL L(500) B (50) 
INTEGER MI(50),MT(50) 
SUM=0. 
DO 1 I=1,N 
NN=1000060 *p (I) 
CALL SIMULA(NN, 
1 SUM=SUM+X 
PRO=SUM/1000000. 
RETURN 
END 


SUBROUTINE FOR SIMULATION 
SUBROUTINE SAN SON’ epitlpeuipe 2 A, LITE) 


Remsen (500) ,R 
INTEGER NE cl Ps 


SMTS MDT +M—1 
DO 11 I=1,NN 


EST=0 


a Ss 
CALL LEXPN(IX,RN,M, 
J J=0 


lee 7 | 0) 


DO 111 J=MIT,MMT 


mm 
= 


JJ=JJI+1 
= ES eerie N 
oLIMe) x 


JJ) /L (J 
GAL 
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BY Beg 


SUBROUTINE CONTRO 
peg lk ee RO) 


LA I(1) ,HT(1) , xX, 1X, TIME) 


SUBROUTINE OU s 
pee oO} 





C 
. Piero Neko lib ols. ot 


WRIT P(1)« (L (2) -2=1, 420) 
sr 


(6 
; 
DO 
UP 
IT 
102 CONTINUE 
WRITE (6, 
RETURN 


c 
CeO sete TATSMENTS 


tHI(T) <1 
eae (bo) -JeNT DO; Mrur) 
TIME, PRO 


6co FORMAT (1 teks) THE SYSTEM IS :') 
601 FORMAT SON KYU BIT) =! P10.7,' * ©, fF EXP L=t,5 
1 B10 81 3X) 03 (//2 £5 (210. 4, 3x) 
6C2 FORMAT('0',9X,'+ ',FPT0.7,' ® 1," EXP L=',5(F1 
10.4, 3X) 65 (7720x ,5(F10.4,3x) )) 
603 FORMAT | O",SX,* PROBABILITY OF SYSTEM SURVIVAL AT 
*tT='°P1d.5,' r§ ', F170. 7) 
END 
INPUT FOR THE PROGRAM : 
TIME 15%, 
#20 Ww 5 
Sy, ay ay. ie 2 
1.87 6 
a4 8 4 i a4 
. 8 
2.XD 7 ; 
Ee a Pee, Ele =) 
25 25 
3. PD g 
gu 5 Ens a4 geal 
25 25 5 3 
4. TH 6 
P| oe, el ae | ra! 
Wy. 
5. 7H 9 
ol a 4 4 <2 
3.2 ae pa 4 
QUEPUT FRON THE PROGRAM : 
THE SYSTEM IS : 
F(T)= 0.2 * BXP OL = 0. 1 SC mOMLIO. 1 Oe 
peomoeeePr | 6 = O12 0.2 #O.3 0.6 0.5 
0.5 0.5 
MUMS ESXS o£ =.0.1 0.5 0.3. 0.1. 0.1 
5 e538 86 osSem 0.3 
fOy lke EXP & = 0.1 meno "Ot One 
foe SX> t= 0,4. O.2 Oc O.1 0.2 
02 EO AD Ss 0 ae 
PROBABILITY OF SYSTEM SURVIVAL AT T= 15.00000 2S 0.911241 
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E.1 Redundant systems with 
Beh s > 
Hee Selaagbility analysis, 


describe a set of components 


sion. A syst2m is redundant if, 


components involved act as back 


Petcidiuederin tion miqht be 
tanz if the failure of any one 
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Z2ts components can fail without 
Thus redundant systems have the 
the failure 
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réviews some notions 


which are found in the 


system is used to 


Organized to perform some nis- 


3n some fashion, some of the 


UOCE OR Orie: COlpenen so. 


that 2 system is not redun- 


of its components causes she 
E2dugdant =f Ons Os MOre Cf 
Causing the systemito East. 
property that they can suf- 


of some of their components 


(ESARY” Ref.1]) 


set cf exoonentially lived ccenpo- 





We have n components, 2ach ind 


wW 
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D 
ry 
fu 
iD 
AD | 
ct 
rh 
if 
O 
3 
cf 
Hg 
wD 
O 
ct 
ee 7 
MD 
if 
in 
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we want =o compute the probability tha; 


wos first. 


P(j.th component fails first)=P(T) <T. ,yi +1#3) 


jee) 


Se 
= | P(t. >T; pyi Ji#IT)=s) Aye? ds 
0 


={ [1 F. (5) ) Jé, (s) as 

1Fj 

oo aves -\.S 
Cece ds 

i#j ; n 
0 n oe n =) 1445 
= A/C) rf CLA TS ds 
izl | 2 i=l 


n 
= ares d. ). 


Peoevedeneracy at zero (Zaro Distributicn): 


Pete o mm Do the name for {he distribution of & Tandon 


variable that is degenerate at zero. 
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